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Abstract 

We use the Mellin-Barnes representation in order to improve the theoretical estimate of mass 
corrections to the width of light pseudoscalar meson decay into a lepton pair, P — > . The full 
resummation of the terms ln(m 2 /A 2 ) (m 2 /A 2 ) n and (m 2 /A 2 ) n to the decay amplitude is performed, 
where mi is the lepton mass and A ~ m p is the characteristic scale of the P — > 7*7* form factor. 
The total effect of mass corrections for the e + e~ channel is negligible and for the /U + ' [i~ channel 
its order is of a few per cent. 

PACS numbers: 13.25. Cq, 12.38.-t, 12.38.Lg 
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I. INTRODUCTION 



Rare decays of mesons serve as a low-energy test of the Standard Model. Accuracy 
of experiments has increased significantly in recent years. Theoretically, one of the main 
limitations comes from the large distance contributions of the strong sector of the Standard 
Model where the perturbative QCD theory does not work. However, in some important 
cases the result can be essentially improved by relating these poorly known contributions 
to other experimentally known processes. The famous example is the calculation of the 
hadronic vacuum polarization contribution to the anomalous magnetic moment of muon 
(g — 2)^ t where the data of the processes e + e~ — > hadrons and r — > hadrons are essential to 
reduce the uncertainty (see for review 

Rasa 

). It turns out that this is also the case for 
the rare decays of light pseudoscalar mesons into a lepton pair 5|. Interest in these processes 
revived after new precise measurement of the decay 7r — > e + e _ by the KTeV collaboration 
6]. The Standard Model prediction disagrees with the KTeV measurement by 3.3cr, with 
the theoretical accuracy exceeding the experimental one. 



P(3) 




FIG. 1: Triangle diagram for the P — > l + l process with the pseudoscalar meson form factor 
P — > 7*7* in the vertex. 



In the lowest order of QED perturbation theory, the photonless decay of the neutral 
meson, P(q) — > l~(p-) +l + {p+), q 2 = M 2 , p\ = m 2 , (M meson mass, m lepton mass) is 
described by the one-loop Feynman amplitude (Fig. [1]) corresponding to the conversion of 
the meson through two virtual photons into a lepton pair. The normalized branching ratio 
is given by [7, 8, 9J 
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where f3 (q 2 ) = a/1 — Amf/q 2 and the reduced amplitude is 
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with the transition form factor Fp 7 * 7 . (— fc 2 , — g 2 ) being normalized as Fp 7 * 7 *(0,0) = 1. 

The imaginary part of the on-shell amplitude A(q 2 = M 2 ) comes from the contribution 
of real photons in the intermediate state and can be found in a model independent way 8] 



Ihu4(M 2 



\n(y(M 2 )), y(q 2 ) 



2/?(M 2 )^ vy ,n u V ! 1 + P(q 2 )' 
A once-subtracted dispersion relation for the amplitude in Eq. ([2]) is written as 1 [l 

ImA(s) 



A(q 2 )=A(q 2 = 0) + ^- r ds- 
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s [s — q^) 

The second term in Eq. (Tj0) takes into account strong q 2 dependence of the amplitude around 
the point q 2 = occurring due to the branch cut coming from the two-photon intermediate 
state. Integrating Eq. (j3J) for q 2 > 4m 2 one arrives at H Q, Q 



Re„4(g 2 )=„4(g 2 = 0)+-^ 



^n 2 (y(. 2 )) + ^ + Li 2 (-y(g 2 )) 



(5) 



where Li2 (z) = — J* (dt/t) In (1 — t) is the dilogarithm function. 

Usually, the subtraction constant in ([5]), containing the nontrivial dynamics of the process, 
is calculated within different models describing the form factor Fp 7 » 7 «(fc 2 , q 2 ) (e.g. 
12]). However, it has recently been shown in [^j that this constant may be expressed in 
terms of the inverse moment of the transition form factor given in symmetric kinematics of 
spacelike photons, G(t) = Fp 7 * 7 * (t, t) , 



A (q 2 = 0) = 3 In 



m 



3 
2 



t 



1 f°° , G (t) 
- + / dt^^- 
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Here, /i is an arbitrary (factorization) scale. One should note that the logarithmic depen- 
dence of the first term on [i is compensated by the scale dependence of the integrals in the 
brackets. 

The accuracy of these calculations is determined by omitted small power corrections of 
the order 0(^,^ln^) and 0(^,^ln^) in the r.h.s. ©, where A < M p is the 



1 In this derivation it is tacitly assumed that the imaginary part of the off-shell amplitude A(q 2 ) has the 
same form as in §5§ with M 2 substituted by q 2 . 
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characteristic scale of the form factor G(t). The aim of this work is to improve the result (J6]) 
for the amplitude A (q 2 = 0) of the P — > / + /~ decay by taking into account all order mass 
corrections ~ ^, ^ In ^. 



II. MELLIN-BARNES INTEGRAL REPRESENTATION 



We evaluate the amplitude A(q 2 ) following the way used in [141 ] . Let us transform the 
integral in (j2J) to the Euclidean metric ko — > %k±. The corresponding integral is convergent 
due to decreasing of Fp 7 * 7 *(fc 2 , (q — k) 2 ) in the Euclidean region. Then use the double Mellin 
transformation for the meson form factor 

2 W a 2 

(2ni 



F Prr (k 2 ,(q-k) 2 ) = 7 ^- I f dz®{z 1 ,z 2 )(^-\ { 



(7) 



Xk-q) 

where A is the characteristic scale for the form factor, dz = dz\dz 2 , the vector o = (a x , cr 2 ) £ 
IR 2 , and $ (zi, z 2 ) is the inverse Mellin transform of the form factor 



00 poo 

dt\ I dt 2 t x t 2 2 Fpj*^* (ti,t 2 ) 
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which has singularities at Re(^) = 0,-1,... . Introducing Feynman parameters in the 
standard way, the denominator part of the integrand in (|2j) can be written as 

(9) 



{k 2 ) Zl+1 [(k - qf] Z2+1 [(p_ - k) 2 + m 2 ] 



r (3 + Zl + z 2 ) 



r(^ + i)r(^ 2 + i) 
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a\a 2 p 2 ). Then the k— loop integral reduces to 
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Combining all factors we get 



(2m) 
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In the general case, to step further we need to perform the third Mellin transformation 
for denominators containing [yj]. Then, considering the process P — ► / + /~ with mass 



hierarchy m << M < A ~ m p we expand the integral obtained over the ratios of the lepton 
and meson masses to the characteristic scale of the meson form factor A by closing the Mellin 
contours in an appropriate manner. However, in the present study we are interested in the 
amplitude at q 2 = 0. In this limit, the Feynman parameter integrals in (fTUj) can be carried 
out in terms of T-functions, and we obtain the following Mellin-Barnes representation for 
A (q 2 = 0) 



a(o) = -^[ dz(ey 



r(z 1 )r(z 2 )r{z 12 )r(i-2z 12 



;-3 + 2zi 2 )$(zi,z 2 ) 



(2iri) Ja+m? ^ T (3 - z V2 ) 

(11) 

with a in the triangle {xi > 0, x 2 > 0, x 2 + X\ < |} chosen so that the integration path 
a + 2IR 2 does not intersect the T-function singularities given by the polar complex lines (see 
illustration in Fig. [2^) 

Li : {zi = -v) , L 2 : {z 2 = -u} , L 3 : {z x + z 2 = -u} , (12) 
L 4 : {l-2( Zl + z 2 ) = -v}, v = 0,1,2,... 

In ( TTT1) we introduce the notation ^ 2 = m 2 /A 2 , z\ 2 = z\ + z 2 and combine the regular 
expression in the squared brackets. 

In further analysis of the integral fTTTj) we use the technique suggested in [15|]. 2 Following 
this line let us associate the vectors in 2-dimensional space with the coefficients of the 
T-function arguments in the numerator and denominator of the integrand in fTTTj) a\ = 
(1, 0) , a 2 — (0, 1) , 03 = (1, 1) , 04 = (—2, —2) , c\ = (—1, —1). Next, define the vector 

and draw through a the straight line /a = {a; G 1R 2 : (A,x) = (A, a)} with the nor- 
mal vector A. The scalar product is introduced as (x, y) = x±yi + x 2 y 2 . The point 
cr divides /a into two rays / + and l~ so that the pair of directions / + and A yields 
the same orientation of 1R 2 as the pair of coordinate axes x\ and x 2 . The half-plane 
7Ta = {x G M 2 : (A,x) < (A, a)} with boundary /a an d the integration half-space Ha = 
7Ta + iM? = {z G C 2 : Re (A, z) < Re (A, cr)} characterize the domain in the space of inte- 
gration variables z in which the integrand is a decreasing function. 



2 When our study was completed we became aware of the results of the work 16| where similar technique 
of the two-dimensional counter integrals is used. 
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FIG. 2: a) Singularities of the integrand in (jilt) . The semi-planes where the arguments of T 
functions produce singularities are depicted by lines Lj with shadowed bands. The point a charac- 
terizing the integration counter is from the triangle {xi > 0, x 2 > 0, x 2 + x\ < ^} b) Rotation of 
/a allows one to read off the degeneration. 

Now we need to define the divisors given by the condition 

D 1 = U {Lj : Lj n r = 0} , D 2 = U {Lj : Lj n l + = 0} . (14) 

The theorem [ijj] states that in the nondegenerated case (A / and all Oj A) the integral 
like (TTTi) is given by the sum 

A(0)= res [Integrand^ (0)] , (15) 

where res [Integrand^. (0)] is the residue with respect to the system of divisors {Di, D 2 }. 

The integral ( |TT1) corresponds to the degenerate case since 03 and CI4 are parallel to A. In 
this case, one has L2 G D\ and L\ G D%. However, L3 and L4 are parallel to /a and cannot 
to be ascribed to any of divisors. To read off the degeneration we slightly rotate /a with 
respect to the point a in clock- wise and anti-clock- wise directions (Fig. [2)3). Now we have 
crossings of L 3 4 with the rotated line and are able to decide to which divisor they should 
be related. One has two cases 

1) D 1 = {L 2 ,L 4 },D 2 = {Li,L 3 }, (16) 

2) D 1 = {L 2 ,L 3 },D 2 = {L 1 ,L i }. 

We are interested in the intersection of divisors, i.e., intersections of all l[ G Di and all 
l}p G D 2 so that these intersections fl belong to the half-space IIa- 

6 



Another important property is the intersection rank, the number of lines that meet at 
each point. If only two lines and meet at each point z r G -DinZ^nllA (rank 1), then 
one has only simple poles. If the intersection rank is more than one, than one may either 
apply the theory of multiple residues or introduce small e— parameters in the arguments 
of T— functions in such way that all poles become simple ones (like in the dimensional 
regularization method). 

We prefer here the second approach, namely, Li,L 2 ,L% in (|12p meet at the same points 
(— a, — (3) where a, (3 = 0,1,... are positive integers. In order to get rid of this kind of 
degeneracy, we add a small parameter e to the argument V (z 2 ) — > T {z 2 + e) in (ITT]) . Now 
we ready to analyze the poles and their residues. Consider first the case 1) in ffTBT) . We have 
two sets of intersections in D\ H D 2 n 7Ta 



L 2 nLi,L 2 nL 3 



(17) 



which may be parametrized as 



L 2 n Zq : 



z 2 + e = -a, 



L 2 n L, : 



z 2 + e = -a, 



z\ = — P, I Z\ = a — P + e, 

Calculating residues we get two contributions to the integral 

A (0) = Aa (0) + A b (0) 



(18) 



(19) 



A (0) = £ 



-1 



a\(3\ 



+/3+£ r {-a - p - e) T (1 + 2 {a + (3 + e)) 



a,/3=0 

(-3-2(a + p + e)) 



T (3 + a + P + e) 



T(-a)T(-P) 



A b (0) = 

a,/3=0 



-1 



/3 r (1 + 2p) 
r(3 + /3) 



-3 - 2/?) 



$ (at — P + e, -a) 
rf-a) 



(20) 



(21) 



The second case in (fl6l) is reduced to the first one because one has single parameter (£) 
integral. If one would be interested in the expansion in inverse powers of £ one needs to 
consider intersections L 4 fl L 1; L 4 fl L 3 instead of (TT71) . 

By using the representation (jSJ) one may show that the corresponding residues are 

$(-a,-/3) , .c+^m ,„ n > $(2;, -a) 



r(-«)r(-/5) 



(-1)^F^(0,0), 



2-l P (0,a) 
- P 7 * 7 * 



(t,0), (22) 
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where (0, 0) denotes the derivatives of an order of a and (3 in the corresponding 

arguments of the form factor. After these substitutions one sum in (1201) and (l2Tj) may be 
performed with the result 

a to) = - 1 <« 2 )" +e r( "'^i:y +£)) < 3 + 2 (» ■ 



n=0 

a (o) = ± (a- r^^g s < 3 + 2 "> f a" 5 ™ 

n! r (3 + n) T (1 - £ + n) J 



(t), 



n=0 

where we again use G(t) = Fp 7 * 7 *(t,t). Now we expand in e and take the limit e — > with 
the total result 



2\ n 



T (1 + 2n) 



?1=0 



n\ r(l + n)T(3 + n) 



(G (n) (0) [2 + (3 + 2n) (ln4£ 2 - 7 



(23) 



(n + 1) + -0 ( n + - ) - i/j (n + 3) 



(3 + 2n) / dtG {n+1) (t)]nt 



Note that the e _1 poles contained in the intermediate steps of calculations are canceled in 
the final expression. To the lowest orders in £ 2 expansion one gets 

A® (0) = - 3 In - - + 3 / dtG {1) (t) In t , (24) 

^(i)( ) = -£ 2 i G (1) (0) ^51n£ 2 + y^ +5 J™ dtGW(t)lnt . (25) 

The leading order expression (j24j) is in accordance with the result obtained in [5| . In the 
general case it is convenient to convert the sum in ( |23i) into the integral form 



-4(0) = ^ f d y\j l —f- i [( ln4 ^ - 7) (2 +2/) + 2 (1 - y)] G (-4^ 2 ) + 



+ (2 + y) 



dt 



\ntG w (t-4y£, 2 ) 



— -t — 3/ — / — t 

G - — - — —G {-Aye-^e) ~ —G (-4^ 2 ) 



(26) 



Finally, let us consider the form factor we are interested in from a physical point of view 

1 



G(t) 



l + t 

For this form factor from (1241) - (126]) one gets the coefficient of logarithmic term as 

1 



.4(0) 



m 4 



l + Qe- VT^W{l + 8e) ln£ 2 + 0(£°) 



(27) 
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or the first terms of expansion 

a(o) = 1(i+ + o - \ (i + %e + o r)) . (28) 

Thus, one can see that in the realistic case for muon, £ 2 = m^/A 2 ~ m 2 ^/m 2 p ~ 0.02 the 
corrections to the leading order coefficients are of an order of 1% and for an electron pair 
they are negligible. 



III. CONCLUSIONS 



The aim of this paper is to clarify the situation with rare decays of pseudoscalar mesons 
to a lepton pair. The situation became more pressing after recent KTeV E799-II experiment 
at Fermilab in which the pion decay into an electron-positron pair was measured using the 
Kl — > 3ir process as a source of tagged neutral pions 6|. The branching ratio was determined 
to be equal to 

B nl-ld -> e+ e~) = (7.49 ± 0.29 ± 0.25) • 10" 8 . (29) 



The standard mode 



7r — > 77 



gives 



prediction based on the use of CLEO data on the transition form factor 



5] 



which is 3.3a below the KTeV result fl29l) . Therefore, it is extremely important to trace 
possible sources of the discrepancy between the experiment and theory. There are a number 
of possibilities: (1) problems with (statistic) experiment procession, (2) inclusion of QED 
radiation corrections by KTeV is wrong, (3) unaccounted mass corrections are important, 
and (4) effects of new physics. At the moment the last possibilities was reinvestigated. In 



171 ]. the contribution of QED radiative corrections to the 7r° — > e + e~ decay, which must be 
taken into account when comparing the theoretical prediction (|30]) with the experimental 
result, fl29|) was revised. Comparing with earlier calculations (lOj, the main progress is in 
the detailed consideration of the 7*7* — > e + e~ subprocess and revealing of dynamics of large 
and small distances. Occasionally, this number agrees well with the earlier prediction based 
on calculations [lOj and, thus, the KTeV analysis of radiative corrections is confirmed. In 
the present paper, we show that the mass corrections are under control and do not resolve 
the problem. So our main conclusion is that the inclusion of radiative and mass corrections 



is unable to reduce the discrepancy between the theoretical prediction for the decay rate 
fl30|) and experimental result fl29l) . The effects of new physics were considered in (3] where 
the excess of experimental data over theory is explained by the contribution of low mass 
(~ 10 MeV) vector bosons appearing in some models of dark matter. Further independent 
experiment at KLOE, NA48, WASAatCOSY, BES III and other facilities will be crucial for 
resolution of the problem. 
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